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In quantum mechanics, the space-fractional Schrodinger equation provides a natural extension of the 
standard Schrodinger equation when the Brownian trajectories in Feynman path integrals are replaced by 
Levy flights. Here an optical realization of the fractional Schrodinger equation, based on transverse light 
dynamics in aspherical optical cavities, is proposed. As an example a laser implementation of the fractional 
quantum harmonic oscillator is presented, in which dual Airy beams can be selectively generated under off-axis 
longitudinal pumping. © 2015 Optical Society of America 
OCIS codes: 140.3410 , 070.2575, 000.1600 


In the past few decades, several theoretical studies 
have speculated about possible extensions of standard 
quantum mechanics. While the implications and use¬ 
fulness of such theories are still a matter of debate, 
some of their ideas and methods have found interest 
beyond their original framework, with potential appli¬ 
cations to science and technology. An example is pro¬ 
vided by VT -symmetric (non-Hermitian) extension of 
quantum mechanics [1], which is attracting a growing 
interest in optics (see, for instance, [2-6] and references 
therein). Space-fractional quantum mechanics (SFQM), 
introduced in 2000 by Laskin [7—9], provides a different 
generalization of the standard quantum mechanics that 
arises when the Brownian trajectories in Feynman path 
integrals are replaced by Levy flights. SFQM provides an 
interesting fractional physical model in quantum physics 
[10]. Most of current studies on SFQM have focused on 
the mathematical aspects of the theory [11-15] and have 
highlighted subtle issues that arise from the non-local na¬ 
ture of the kinetic operator [13], however very few studies 
have considered possible physical realizations or applica¬ 
tions of SFQM [10,16]. In Ref. [16] a condensed-matter 
realization of SFQM based on Levy crystals was sug¬ 
gested, however its experimental implementation seems 
very challenging owing to the need to specially tailoring 
non-nearest neighbor hopping in the lattice. 

In this Letter we suggest a different route toward 
a physical implementation of SFQM, which is based 
on transverse light dynamics in aspherical optical res¬ 
onators. As an example, we discuss the design of an op¬ 
tical resonator that realizes the fractional quantum har¬ 
monic oscillator [7,9] (also known as massless relativistic 
quantum oscillator [17,18]). Interestingly, we show that 
dual Airy beams [19] can be generated from such an op¬ 
tical cavity. 

The space-fractional Schrodinger equation (SFSE) was 
derived by Laskin is a series of papers using a Feynman 
path integral approach over Levy trajectories [7-9]. In 


the one spatial dimension and assuming % = 1 it reads 

where a is the Levy index (1 < a < 2), V a is a scale 
constant, V(x) is the external potential, and if = if(x, t) 
is the particle wave function. In Eq.(l) the kinetic term 
in the Hamiltonian is represented by the quantum Riesz 
derivative (fractional Laplacian) of order a , which is de¬ 
fined by 
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The ordinary Schrodinger equation is obtained in the 
limiting case a = 2. In Ref. [16], a realization of the 
SFQM based on lattice dynamics in so-called Levy crys¬ 
tals was suggested. Here we propose an optical realiza¬ 
tion of the SFSE based on transverse light dynamics in 
an optical cavity, in which we exploit the properties of 
Fourier optics to realize the fractional Laplacian [20,21]. 
We consider an effective one-dimensional optical res¬ 
onator problem with transverse x spatial direction. Such 
an analysis holds for e.g. an optical cavity with astig¬ 
matic optical elements in which separation of variables 
in the two transverse spatial dimensions x and y can 
be applied and the lowest-order transverse mode in the 
y direction is excited. A schematic of the effective one¬ 
dimensional Fabry-Perot optical resonator is depicted in 
Fig. 1(a). It comprises two flat end mirrors and two con¬ 
verging lens of focal length / in a typical 4/ configura¬ 
tion, with two thin phase masks with transmission func¬ 
tions ti(x) = exp[—if(x)/2] and ^(x) = exp[—iV(x)/2], 
where 

f( x ) = p\x\ a . (3) 

In the scalar and paraxial approximations, wave propa¬ 
gation at successive transits inside the optical cavity can 
be readily obtained from the generalized Huygens inte¬ 
gral by standard methods [22]. Neglecting finite aperture 
effects and indicating by F^ n \x) = i n ip( n \x) the field 
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exp[— iV(x) — if(£)] entering in the kernel K is ex¬ 
panded as exp [—iV(x)—if(£)\ ~ 1 — iV{x) — if{£), using 
Eqs.(2,3) one readily obtains 


V> (n+ 1 ) (x) 
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where we have set V a = /3(A//27r) a . If we introduce 
the time variable £, normalized to the round-trip transit 
time Tr of photons in the cavity, after setting ^(x, t) = 
i/;( n=t )(x), (d'lp/dt) ~ ^ n+1 \x) —^ n \x), the temporal 


evolution of the transverse field envelope if(x,t) at plane 
7 is thus governed by the SFSE 
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Note that at the Fourier plane q the Fourier transform 
field 
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Fig. 1. (Color online) (a) Schematic of an optical resonator 
in a 4-/ configuration that realizes the SFSE. (b) Schematic 
of the optical resonator that realizes the fractional quantum 
harmonic oscillator. The phase mask t 2 [x) and flat end mir¬ 
ror in (a) are replaced by a spherical mirror. The pumping 
scheme to selectively excite higher-order dual-Airy TEM n 
modes is shown in the lower panel. 


of ip(x,t) is observed, which satisfies the SFSE in ’mo¬ 
mentum space’ [ 11 ] 
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envelope of the progressive wave at the reference plane 7 
in the cavity and at the n-th round-trip, in the absence 
of gain/loss elements one can write 

\x) = f d9K(x, 0)^(0), (4) 


where the kernel K of the integral transformation is 
given by 


K(x,6) 
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Therefore the transverse modes and resonance frequen¬ 
cies of the resonator of Fig. 1(a) correspond to the eigen¬ 
functions and corresponding energies of the SFSE with 
Levy index a in the external potential V{x). 

As an example, let us consider the optical realization of 
the fractional quantum harmonic oscillator [7,9], which 
corresponds to a parabolic external potential V (x) oc x 2 . 
This case can be readily implemented by substitution of 
the phase mask and flat mirror at plane 7 by a simple 
spherical mirror of radius of curvature R [see Fig.1(b)], 
corresponding to the potential 

On r 

V(x) = —x 2 . (10) 


and A = 2iv/k = cjv is the optical wavelength. The 
transverse modes TEM n and corresponding resonance 
frequencies of the resonator are obtained from the eigen¬ 
functions 'ipn(x) and corresponding eigenvalues a n of the 
Fredholm equation a n 'ip n (x) = f d0K(x,0)if; n (6) [ 22 ]. 
Note that, in the absence of the phase masks, i.e. for 
V(x) = f(x) = 0 , the 4/ resonator is a self-imaging 
cavity, i.e. K{x,6) = S(x — 0), so that any field dis¬ 
tribution is self-reproduced at each round trip. Here 
we are interested in the limiting case V(x),f(x) -T 0 , 
so that the field is nearly self-imaged and undergoes 
slow changes at each round trip. If the exponential term 


The transverse modes (eigenfunctions) and correspond¬ 
ing frequencies of the resonator can be at best deter¬ 
mined from the SFSE eigenvalue equation in momentum 
space [Eq.(9)], i.e. at the Fourier plane q. In fact, after 
setting 4>(x,t) = cj)(x) exp(— iEt), from Eqs.(9) and (10) 
one obtains the eigenvalue equation 
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A particularly interesting case which admits of an an¬ 
alytical solution is the one corresponding to a = 1 , 
which is realized by using a one-dimensional refractive 
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axicon in the Fourier plane q . This limiting case of frac¬ 
tional quantum oscillator realizes the so-called massless 
relativistic harmonic oscillator, which is studied e.g. in 
Refs. [17,18]. For a = 1 Eq.(ll) reduces to the Airy 
equation, with eigenfunctions (transverse modes TEM n ) 
and corresponding energies E n given by 

4>n(x) = Ai (k|x| + r n ) (12) 

E - = (13) 

where Ai(x) is the Airy function, k = (2tt/ 3R/Xf 2 )^ 1 / 3 ^ 
and where r n are the zeros of ( dAi/dx ) (for n even) or 
of Ai(x) (for n odd). Note that, since r n < 0 , one has 
E n > 0. Interestingly, for large enough mode number n 
the transverse modes TEM n , as described by Eq.(12), 
mainly correspond to two displaced Airy beams that ac¬ 
celerate into opposite directions, the so-called dual Airy 
beams [19]. In an optical experiment, such higher-order 
modes can be selectively excited in a laser device un¬ 
der off-axial longitudinal pumping [23]. To this aim, let 
us introduce a thin gain medium close to the spheri¬ 
cal mirror, see Fig. 1(b). We assume that the medium is 
longitudinally-pumped by two off-axis Gaussian beams 
of spot size w p and transversely displaced by 2 p. Taking 
into account the presence of the gain medium and the 
output coupling of the cavity, the round-trip propagation 
of the intracavity field ^ n \x ) at plane 7 is modified as 
follows [compare with Eq.(4)] 

^ nJrl \x) = V^Texp^x)] J d0K(x 1 0)'ip^ n \0) (14) 

where T is the transmittance of the output coupler, g(x) 
is the transversely-dependent round-trip gain in the ac¬ 
tive medium, and the kernel K is given by Eq.(5). The 
gain parameter g(x) is given by g(x) = 2crN(x)l, where 
l is the thickness of the gain medium, a is the stimu¬ 
lated cross section of the lasing transition and N(x) is 
the population inversion in the medium. For the off-axis 
longitudinal pumping scheme shown in Fig. 1(b) and ne¬ 
glecting saturation of the gain, we may assume 


by the Fox-Li iterative method for parameter values 
/ = 1 cm, A = 1064 nm (Nd:YAG laser), T = 97%, 
R = 50 cm, f3 = 5 x 10 -3 /rm -1 , w p = \[2 x 20 pm 
and for increasing values of p. For p = 0, i.e. for on-axis 
pumping [Fig.2(a)], the lasing mode is the lowest-order 
transverse mode TEMo, corresponding to ro — —1.05 
in Eq.(12). At p = 170 pm [Fig.2(b)], the lasing mode 
is the transverse mode TEM6 with mode index n = 6, 
corresponding to vq ~ —6.16. Finally, at p = 256 pm 
[Fig.2(c)] the lowest-threshold mode is the transverse 
mode TEM 16 with mode index n = 16, correspond¬ 
ing to ri6 — —11.48. The threshold values of lasing in 
three cases are goth — 0.028 in Fig.2(a), goth — 0.075 
in Fig.2(b) and goth — 0.098 in Fig.2(c). Therefore the 
laser system schematically shown in Fig. 1(b) provides an 
optical realization of the fractional quantum harmonic 
oscillator, in which the different modes of the oscillator 
can be selectively excited by suitable off-axis pumping. 
For large off-axis pumping, the lasing modes correspond 
to dual Airy beams, and thus our system provides an 
interesting example (additional to those previously re¬ 
ported in Refs. [23,25]) where Airy-type beams can be 
directly generated from a laser oscillator. It should be 
finally noted that the above analysis does not consider 
mode confinement in the other transverse y direction, 
however such a confinement can be readily obtained us¬ 
ing an astigmatic optical cavity. Namely, let us consider 
the optical cavity of Fig. 1(b) in which the two converg¬ 
ing lenses (as well as the phase mask at the Fourier plane 
q) are astigmatic and acts along the spatial x direction 
solely, whereas the end spherical mirror is circularly sym¬ 
metric and so it acts in both x and y directions. In 
this case the resonator TEM n?m transverse modes can 
be readily obtained by separation of variables, i.e. they 
factorize as the products of the modes of two effective 
one-dimension resonators in the x and y directions. For 
example, at the Fourier plane q the transverse modes 
TEM„, m are given by 4> n ,m(x,y) = 4> n (x)Y m (y), where 
<p n (x) are the dual Airy modes given by Eq.(12) whereas 

Y m (y) = H m exp (-A “ ( 16 ) 
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where w p is the pump beam spot size and go is the 
peak gain parameter. The lowest-thresold lasing mode 
and corresponding laser threshold go = goth can be nu¬ 
merically computed from Eq.(14) by a standard Fox- 
Li iterative method [24]. As the displacement 2 p be¬ 
tween the two pumping beams is increased, higher- 
order TEM n dual-Airy transverse modes of the res¬ 
onator with increasing mode number n can be selec¬ 
tively excited. This is clearly shown in Fig.2. The figure 
shows the numerically-computed lowest-threshold trans¬ 
verse modes, both in the near-field plane 7 (left panels) 
and in the Fourier plane q (right panels), as obtained 


are standard Gauss-Hermite modes of index m with 
beam spot size w y and radius of curvature R y . The values 
of w y and R y at the Fourier plane q can be readily com¬ 
puted by standard Gaussian beam analysis of the res¬ 
onator in the ^/direction and read 


- - m- o”, 

/ 4R — 12 / 

(17) 

/ £-4/ 

Ry = 2 R- 6 /. 


(18) 


For a two-beam longitudinal pumping with 2 p off- 
axis displacement in the x direction solely, i.e. for 
a two-dimensional optical gain distribution g(x,y ) = 
go exp[— 2 (x - p) 2 /w 2 p - 2y 2 /w 2 p } + +00 exp[- 2 (a; + 
p) 2 /w 2 — 2 y 2 /w 2 \, the lowest order transverse mode 
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m = 0 in the y direction is excited. In this case, as p 
is increased from zero, higher order dual-Airy transverse 
modes TEM n? o in the x direction are excited, and the re¬ 
sulting two-dimensional intensity distributions of lasing 
modes are those depicted in Fig.3. 


near-field plane y Fourier plane q 



Fig. 2. (Color online) Numerically-computed intensity dis¬ 
tribution of the lowest-threshold lasing modes in the optical 
cavity of Fig. 1(b) for increasing values of the off-axis pump 
displacement p. The left panels show to the transverse inten¬ 
sity distributions in the near-field plane 7 , whereas the right 
panels are the corresponding intensity distributions in the 
Fourier plane q. (a) p — 0 (on-axis pumping), (b) p = 170 pm, 
and (c) p — 256 pm. The other parameter values are given in 
the text. The dashed curves in the right panels, almost over¬ 
lapped with the solid ones, show the intensity distributions 
of the transverse modes TEM n , as predicted by Eq.(12), for 
n = 0 [in (a)], n — 6 [in (b)], and n— 16 [in (c)]. 



Fig. 3. (Color online) Two-dimensional intensity distribu¬ 
tions of the lasing modes at the Fourier plane q in the astig¬ 
matic cavity of Fig. 1(b), taking into account confinement in 
the other spatial direction y provided by the spherical end 
mirror. For increasing values of the the pump displacement p 
higher-order TEM n o transverse modes are excited: (a) p — 0 
(on-axis pumping), (b) p — 170 pm, and (c) p — 256 pm. 

In conclusion, an optical implementation of the frac¬ 
tional Schrodinger equation has been proposed, which 
is based on transverse laser dynamics in aspherical op¬ 
tical resonators. As an example, a laser design for the 


generation of dual-Airy states of a fractional quantum 
harmonic oscillator has been presented. The ability in 
optics to implement fractional pseudo-derivative opera¬ 
tors could be further exploited to simulate other quan¬ 
tum models involving fractional operators, such as the 
relativistic Salpeter wave equation [26,27]. Our results 
indicate that optics can provide a laboratory tool where 
fractional models developed in quantum physics can be¬ 
come experimentally accessible. Fractional models could 
also provide new interesting ways to mold diffraction of 
light and to generate novel beam solutions [28]. 

References 

1. C.M. Bender, Rep. Prog. Phys. 70, 947 (2007). 

2. R. El-Ganainy, K. G. Makris, D. N. Christodoulides, and 
Z.H. Musslimani, Opt. Lett. 32, 2632 (2007). 

3. K. G. Makris, R. El-Ganainy, D. N. Christodoulides, and 
Z. H. Musslimani, Phys. Rev. Lett. 100, 103904 (2008). 

4. C.E. Riiter, K.G. Makris, R. El-Ganainy, D.N. 
Christodoulides, M. Segev, and D. Kip, Nature Phys. 
6 , 192 (2010). 

5. S. Longhi, Phys. Rev. A 82, 031801 (2010). 

6 . A. Regensburger, C. Bersch, M.-A. Miri, G. On- 
ishchukov, D.N. Christodoulides, and U. Peschel, Nature 
488, 167 (2012). 

7. N. Laskin, Phys. Lett. A 268, 298 (2000). 

8 . N. Laskin, Phys. Rev. E 62, 3135 (2000). 

9. N. Laskin, Phys. Rev. E 66 , 056108 (2002). 

10. R. Herrmann, Fractional Calculus: An Introduction for 
Physicists (World Scientific, 2011), chaps. 9 and 10. 

11. J. Dong and M. Xu, J. Math. Phys. 48, 072105 (2007). 

12. E. Capelas de Oliveira and J. Vaz, J. Phys. A 44,185303 

( 2011 ). 

13. Y. Luchko, J. Math. Phys. 54, 012111 (2013). 

14. M. Zaba and P. Garbaczewski, J. Math. Phys. 55, 
092103 (2014). 

15. J.D. Tare, J. Perico, and H. Esguerra, Physica A 407, 
43 (2014). 

16. B. A. Stickler, Phys. Rev. E 88 , 012120 (2013). 

17. K. Kowalski and J. Rembielinski, Phys. Rev. A 81, 
012118 ( 2010 ). 

18. J. Lorinczi and J. Malecki, J. Diff. Eq. 253, 2846 (2012). 

19. C.-Y. Hwang, D. Choi, K.-Y. Kim, and B. Lee, Opt. 
Expr. 18, 23504 (2010). 

20. H. Kasprzak, Appl. Opt. 21, 3287 (1982). 

21. J.A. Davis, D.A. Smith, D.E. McNamara, D.M. Cottrell, 
and J. Campos, Appl. Opt. 40, 5943 (2001). 

22. A.E. Siegman, Lasers (University Science, Mill Valley, 
Calif., 1986), Chaps. 15, 16 and 20. 

23. S. Longhi, Opt. Lett. 36, 716 (2011). 

24. A. O. Fox and T. Li, Bell Syst. Tech. J. 40, 453 (1961). 

25. G. Porat, I. Dolev, O. Bariev, and A. Arie, Opt. Lett. 
36, 4119 (2011). 

26. E.E. Salpeter, Phys. Rev. 87, 328 (1952). 

27. P. Cea, P. Colangelo, G. Nardulli, G. Paiano, and G. 
Preparata, Phys. Rev. D 26, 1157 (1082). 

28. J.C. Gutierrez-Vega, Opt. Lett. 32, 1521 (2007). 


4 




































References with full titles 


1. C.M. Bender, Making sense of non-Hermitian Hamiltonians, 
Rep. Prog. Phys. TO, 947 (2007). 

2. R. El-Ganainy, K. G. Makris, D. N. Christodoulides, and Z.H. 
Musslimani, Theory of coupled optical VT symmetric structures, 
Opt. Lett. 32 , 2632 (2007). 

3. K. G. Makris, R. El-Ganainy, D. N. Christodoulides, and Z. H. 
Musslimani, Beam Dynamics in VT-Symmetric Optical Lattices, 
Phys. Rev. Lett. 100 , 103904 (2008). 

4. C. E. Riiter, K.G. Makris, R. El-Ganainy, D.N. Christodoulides, 
M. Segev, and D. Kip, Observation of paritytime symmetry in op¬ 
tics, Nature Phys. 6, 192 (2010). 

5. S. Longhi, VT symmetric laser-absorber, Phys. Rev. A 82 , 
031801 (2010). 

6. A. Regensburger, C. Bersch, M.-A. Miri, G. Onishchukov, D.N. 
Christodoulides, and U. Peschel, Parity-time synthetic photonic 
lattices, Nature 488, 167 (2012). 

7. N. Laskin, Fractional Quantum Mechanics and Levy Path Inte¬ 
grals, Phys. Lett. A 268 , 298 (2000). 

8. N. Laskin, Fractional Quantum Mechanics, Phys. Rev. E 62 , 
3135 (2000). 

9. N. Laskin, Fractional Schrodinger equation, Phys. Rev. E 66, 
056108 (2002). 

10. R. Herrmann, Fractional Calculus: An Introduction for Physi¬ 
cists (World Scientific, 2011), chaps. 9 and 10. 11. J. Dong and 
M. Xu, Some solutions to the space fractional Schrodinger equa¬ 
tion using momentum representation method, J. Math. Phys. 48, 
072105 (2007). 

12. E. Capelas de Oliveira and J. Vaz, Tunneling in fractional 
quantum mechanics, J. Phys. A 44,185303 (2011). 

13. Y. Luchko, Fractional Schrodinger equation for a particle mov¬ 
ing in a potential well, J. Math. Phys. 54, 012111 (2013). 

14. M. Zaba and P. Garbaczewski, Solving fractional Schrodinger- 
type spectral problems: Cauchy oscillator and Cauchy well, J. 
Math. Phys. 55, 092103 (2014). 

15. J.D. Tare, J. Perico, and H. Esguerra, Transmission through 
locally periodic potentials in space-fractional quantum mechanics, 
Physica A 407 , 43 (2014). 

16. B. A. Stickler, Potential condensed-matter realization of space- 
fractional quantum mechanics: The one-dimensional Levy crystal, 
Phys. Rev. E 88, 012120 (2013). 

17. K. Kowalski and J. Rembielinski, Relativistic massless har¬ 
monic oscillator, Phys. Rev. A 81 , 012118 (2010). 

18. J. Lorinczi and J. Malecki, ” Spectral properties of the massless 
relativistic harmonic oscillator”, J. Diff. Eq. 253, 2846 (2012). 

19. C.-Y. Hwang, D. Choi, K.-Y. Kim, and B. Lee, Dual Airy 
beams, Opt. Expr. 18 , 23504 (2010). 

20. H. Kasprzak, Differentiation of a noninteger order and its op¬ 
tical implementation, Appl. Opt. 21 , 3287 (1982). 

21. J.A. Davis, D.A. Smith, D.E. McNamara, D.M. Cottrell, and 
J. Campos, Fractional derivatives-analysis and experimental im¬ 
plementation, Appl. Opt. 40 , 5943 (2001). 

22. A.E. Siegman, Lasers (University Science, Mill Valley, Calif., 
1986), Chaps. 15, 16 and 20. 

23. S. Longhi, Airy beams from a microchip laser, Opt. Lett. 36, 
716 (2011). 

24. A. O. Fox and T. Li, Resonant modes in a maser interferom¬ 
eter, Bell Syst. Tech. J. 40 , 453 (1961). 

25. G. Porat, I. Dolev, O. Bariev, and A. Arie, Airy beam laser, 
Opt. Lett. 36, 4119 (2011). 

26. E.E. Salpeter, Mass Corrections to the Fine Structure of 
Hydrogen-Like Atoms, Phys. Rev. 87 , 328 (1952). 

27. P. Cea, P. Colangelo, G. Nardulli, G. Paiano, and G. Preparata, 
WKB approach to the Schrodinger equation with relativistic kine¬ 
matics, Phys. Rev. D 26, 1157 (1082). 

28. J.C. Gutierrez-Vega, Fractionalization of optical beams: I. Pla¬ 
nar analysis, Opt. Lett. 32 , 1521 (2007). 


5 


